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Abstract 

We study the asymptotic behavior of the emptiness formation probability for large spin 
strings in a translation invariant quasifree nonequilibrium steady state of the isotropic 
XY chain. Besides the overall exponential decay, we prove that, out of equilibrium, the 
exponent of the subleading power law contribution to the asymptotics is nonvanishing 
and strictly positive due to the singularities in the density of the steady state. 
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1 Introduction 

In this note, I propose to enlarge upon the study started in Aschbacher [6] of the asymptotic 
behavior of a special and important correlator, the so-called emptiness formation probability 
(EFP). Written down in the framework of a spin system over the two-sided discrete line, the 
EFP observable is given by 

n 

A n = ]Jpi, (1) 

i=l 
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where pi denotes the projection p := (1 - cr 3 )/2 at site i onto the spin down configuration of 
the spin a 3 in the 3-direction. For a given state u of the spin system, the probability that a 
ferromagnetic string of length n is formed in this state is thus expressed by 



Due to the existence of the Jordan-Wigner transformation which maps, in a certain sense, 
spins onto fermions, the EFP has been heavily studied for states of the XY chain whose 
formal Hamiltonian is given in Remark [3] below. As a matter of fact, this model becomes a 
gas of independent fermions under the Jordan-Wigner transformation, and it is thus ideally 



The large n behavior of the EFP in the XY chain has already been analyzed for the cases 
where the state u is a ground state or a thermal equilibrium state at positive temperature. In 
both cases, the EFP can be written as the determinant of the section of a Toeplitz operator 
with scalar symbol. Since the higher order asymptotics of a Toeplitz determinant is highly 
sensitive to the regularity of the symbol of the Toeplitz operator, the asymptotic behavior of 
the ground state EFP is qualitatively different in the so-called critical and noncritical regimes 
corresponding to certain values of the anisotropy and the exterior magnetic field of the XY 
chainj^] It has been found in Shiroishi et al. [TS1 that the EFP decays like a Gaussian in 
one of the critical regimesF] In a second critical regime and in all noncritical regimes, the 
EFP decays exponentiallyJjThese results have been derived by using well-known theorems 
of Szego, Widom, and Fisher-Hartwig, and the yet unproven Basor-Tracy conjecture and 
some of its extensions, see Widom [19] and Bottcher and Silbermann [TTJ |T2). On the 
other hand, in thermal equilibrium at positive temperature, the EFP can again be shown to 
decay exponentially by using a theorem of Szego, see for example Shiroishi et al. [T6) and 
Franchini and Abanov 1131. 



In contrast, for the case where u is the nonequilibrium steady state (NESS) constructed 
in Aschbacher and Pillet [9]0 the EFP can still be written as a Toeplitz determinant, but 
now, the symbol is, in general, no longer scalar. Due to the lack of control of higher order 
determinant asymptotics in Toeplitz theory with nontrivial irregular block symbols, I started 
off by studying bounds on the leading asymptotic order for a class of general block Toeplitz 
determinants in Aschbacher [6j. It turned out that suitable basic spectral information on the 
density of the state is sufficient to derive a bound on the rate of the exponential decay of 
the EFP in general translation invariant fermionic quasifree states. This bound proved to 
be exact not only for the decay rates of the EFP in the ground states and the equilibrium 
states at positive temperature treated in Abanov and Franchini p] [T3) and Shiroishi et al. 

1 Low-dimensional magnetic systems are also heavily studied experimentally, see for example Sologubenko 
etal. 02. 

2 l.e., in (T5) below, the parameters 7 and A, respectively. 

3 With some additional explicit numerical prefactor and some power law prefactor, see Shiroishi et al. (T6] 
and references therein. 

4 ln contrast to the noncritical regime, there is an additional power law prefactor in the second critical regime 
whose exponent differs from the one in the first critical regime, see Abanov and Franchini [T|[T3]. 
5 And in Araki and Ho [5] for 7 = A = 0. 



P(n)=u(A n ). 



(2) 



suited for rigorous analysis^ 
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[16], but it will also do so for the nonequilibrium situation treated here exhibiting the so- 
called left mover-right mover structure already found in Aschbacher [7] and Aschbacher and 
Barbaroux [8j for nonequilibrium expectations of different correlation observablesj^] Hence, 
given this exponential decay in leading order which parallels qualitatively the behavior in 
thermal equilibrium at positive temperature, one may wonder whether there is some char- 
acteristic signature of the nonequilibrium left at some lower order of the EFP asymptotics. It 
turns out, and this is the main result of this note, that, in contradistinction to the leading order 
contribution, the subleading power law contribution to the large n asymptotics of the EFP in 
Fisher-Hartwig theory is sensitive to the singularity of the symbol of the underlying Toeplitz 
operator, and it has a strictly positive exponent if and only if the system is truly out of equi- 
librium. This may be interpreted as the manifestation, in subleading order, of the long-range 
nature of the underlying formal effective Hamiltonian of the NESS^This connection is not 
made precise here, though, but it is left to be studied in greater detail elsewhere. 

Section [2] contains the setting and Section |3]the main assertion. The reader not familiar 
with quasifree states on CAR algebras and/or with Toeplitz theory may consult the Appendix 
where definitions and basic facts are collected. 



2 Nonequilibrium setting 

In this section, I shortly summarize the setting for the system out of equilibrium used in 
Aschbacher and Pillet [9j. In contradistinction to the presentation there, I will skip the formu- 
lation of the two-sided XY chain as a spin system and rather focus directly on the underlying 
C*-dynamical system structure in terms of Bogoliubov automorphisms on a selfdual CAR 
algebra as introduced by Araki [4\^\ 

For some iV e {0} u N, the nonequilibrium configuration is set up by cutting the finite 
piece 

Z s := {x G Z | - N < x < N} (3) 

out of the two-sided discrete line Z. This piece will play the role of the confined sample, 
whereas the remaining parts, 

Z L := {x G Z\x < -(N + 1)}, (4) 

Z R := {x G Z\x > N + 1}, (5) 

will act as infinitely extended thermal reservoirs at different temperatures to which the sample 
will be suitably coupled. 

6 This has already been noted in Aschbacher [6J. 

7 See Remark 3 in Aschbacher and Pillet [9]. This effective Hamiltonian is to be understood on a formal level 
only. It has been shown by Matsui and Ogata [1 4] that there exists no dynamics on the Pauli spin algebra w.r.t. 
which this NESS is a KMS state. 

8 For an introduction to the algebraic approach to open quantum systems, see also for example Aschbacher 

etal. ma. 
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We first specify the observables contained in the system to be considered. 

Definition 1 (Observables) Let 5(f)) denote the fermionic Fock space built over the one- 
particle Hilbert space of wave functions on the discrete line, 

h:=£ 2 (Z). (6) 

With the help of the usual creation and annihilation operators a*(f), a(f) e £(#(())) for any 
f e rjFlfr/e complex linear mapping B : f)® 2 -> £(#(*))) is defined, for F ■= [fi, f 2 } e f)® 2 , by 



B(F):=a*(f 1 ) + a(f 2 ). (7) 
Moreover, the antiunitary involution J : f)® 2 -> f)® 2 /'s g/Ven, for a// /i,/ 2 e fj, by 

J[/i,/ a ]:=[/ a ,/i]- (8) 

7r;e observables are described by the selfdual CAR algebra over (j® 2 w/Yr> involution J 
generated by the operators B(F) e £(#(f))) for all F e f)® 2 . / denote this algebra by 

2t:=2t([)® 2 , J)F| 



The time evolution is generated as follows. 

Definition 2 (Dynamics) Let A e K, and let u e £(fj) foe tf?e translation operator defined by 
(uf)(x) := f(x - 1) for all f e f) and all x e Z. 7fre coupled and the decoupled one-particle 
Hamiltonians h, h g £(fj) are defined by 

h := Re(w) + A, (9) 
fr : = / l _( t , L+t;jl ) > (10) 

respectively where the decoupling operators v L ,v R e £°(f)) have £r/e form 

UL := Re^-^W"), (11) 
Vfl := Re(VW (JV+1) ), (12) 

anc/ tf?e projection p e £°(rj) /'s g/Ven by p f := {S ,f)5 for all f e rjj^] For a// 4 g E, 
£r/e coupled and the decoupled time evolutions are the Bogoliubov * -automorphisms r*, r * g 
Aut(2l) defined on the generators B(F) e 21 w/f/7 F g t)® 2 by 

t\B(F)) := B(^ h(S ' h) F), (13) 
r *(F(F)) := 5(e i ' (?i "®- /lo) F). (14) 



9 The bounded operators on the Hilbert space H are denoted by C{%). 

10 The concept of a selfdual CAR algebra has been introduced and developed by Araki [2 3J. Here, it is just 
a convenient way of working with the linear combination {Jj. 

11 1 write £°(H) for the finite rank operators on the Hilbert space H. Moreover, S x e t) for x e Z denotes the 
Kronecker function. Finally, for an operator A e C(H), the real part is Re(A) ■= (A + A*)/2. 
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Remark 3 As mentioned above, this model has its origin in the XY spin chain whose formal 
Hamiltonian is given by 



^=-iE{( i +7)-f ) -r i) +(i-7) 

xez 



(15) 



where 7 e (-1, 1) denotes the anisotropy, A e R the external magnetic field, and the Pauli 
basis of C 2x2 reads 



<TQ 
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" 1 





1 




1 







i 




-1 



(16) 



The Hamiltonian h from (|9) corresponds to the isotropic XY chain, i.e. to the case where 

7 = 0. 

The left and right reservoirs carry the inverse temperatures (3 L and (3 R , respectively. Pour 
fixer les idees, we assume w.l.o.g. that they satisfy 



< p L < Pr < OO. 



(17) 



We next specify the state in which the system is prepared initially. It consists of a KMS state 
at the corresponding temperature for each reservoir, and, w.l.o.g., of the chaotic state for the 
sample. For the definition of quasifree states, see Appendix [A} 

Definition 4 (Initial state) The initial state u e is the quasifree state specified by the 
density S e £(f)® 2 ) of the form 



So '■— s o,- © s o,+) 
where the operators s ,± e £(rj) sre defined by 

so, ± := (1 + e^)- \ 
and k e £(fj ~ f) L © fj 5 © fj fl ) is given by 

k ■= /3 L h L © © /3 R h R . 



(18) 



(19) 



(20) 



Here, for a = L,S,R, I used the definitions t) a := £ 2 (Z a ) and h a ■= i* a hi a , where i a : rj a ->■ fj 
is the natural injection defined, for any f e fj a , by i a ({f(y)} ye % a )(x) := f(x) if x e Z w , and 
zero otherwise. 

The following definition is due to Ruelle [T5]. 

Definition 5 (NESS) /\ A/ESS associated with the C* -dynamical system (21, r) having the 
initial state u e £(21) /'s a wea/<-* limit point forT -> 00 o/' ne^ 



T > 



(21) 
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In the model specified by the Definitions [T|[2j and[4} we get the following NESS using the 
scattering approach of Ruelle [15]. 

Theorem 6 (XY NESS) There exists a unique quasifree NESS u e Q(2l) w.r.t. the initial 
state u e (2(21) and the coupled dynamics r* g Aut(2l) whose density S e C{H) reads 

S = s.® s+, (22) 

where the operators s± e £($) act in momentum space rj := L 2 (T) as multiplication by 

s ± (e ifc ):=i(l±0 ± (e ife )), (23) 

and the functions g± : T ->■ (-1, 1) are defined by 

g±(e ik ) :=tanh[i(^±sign(sinA;)(5)(A + cosife)]. (24) 

Here, we set (3 := + L )/2 and 5 := (/3 fl - /3 L )/2, and the sign function sign : R -> {±1} 
is defined by sign(x) := 1 /fx > 0, anc/ sign(x) := -1 if x < 0. 

Proof. See Aschbacher and Pillet [9]. □ 



The main object of our study is the following. 
Definition 7 (NESS EFP) Letn e N. The EFP observable A n e 21 is defined by 



2n 



A n :=\[B{F i ), (25) 



i=i 



where, for alii e N, the form factors F t e f)® 2 are given by 

F 2i -i := u i @u i G 1 , (26) 
F 2i := u l ®u l G 2 , (27) 

and the initial form factors d, G 2 e f)® 2 look like 

Gl :=JG 2 :=[0,5 ]. (28) 

Moreover, the expectation value P : N ->■ [0, 1] of the EFP observable A n e 21 in the NESS 
u> g £(21) /'s denoted by 



12 



P(7i):=w(40. (29) 



12 As for the name EFP, note that A n = f]" =1 a^a*, and that, for £?„ := ]T" =1 a,, we have 

71 

< P(n) - w(B n B;) < ||B„|| 2 < Y[ H^ll 2 = 1. 

i=l 
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The next assertion states the main structural property of the EFP correlation function. 



For the basic facts of Toeplitz theory, see Appendix p2. 



Proposition 8 (EFP determinantal structure) The NESS EFP P : N -> [0, 1] is given by 
the determinant of the finite section of the Toeplitz operator T[s_] e C(£ 2 (N)), 



P(n) = det(T n [s_]). 



(30) 



Proof. Proceeding as in Aschbacher and Barbaroux [8], we have that, on one hand, the 
skew-symmetric EFP correlation matrix Q n e iR 2nx2ri , defined, for i,j = 1, 2n, by 



n,ij 



'u{B(Fi)B{Fi)), \U<j, 
o, if % = j, 

- u(B(Fj)B(Fi)), \U>j, 



(31) 



where F t e P) 02 for % e N are the form factors from Definition |7J relates to the EFP as 

P(n) = pf(fl n ), 
and, on the other hand, that it has the Toeplitz structure 



(32) 



(33) 



Here, a P e L^ x2 (T) is the block symbol of the Toeplitz operator T[a P ] e C(£l(N)) which I 
computed in Aschbacher [6) to be of the form a P = (s - p)a x , where s = s_ © s + is the 
density of the NESS u e Q(2l) in momentum space and p = (1 - cr 3 )/2. Theorem [6] then 
implies that, in the present case, the symbol has the form 



a-p 



s_ 
s+ - 1 



Hence, there exists anfi e 0(2n) with det(i?) = (-l)"^- 1 )/ 2 s.t., using Lemma p9 
reduce the block Toeplitz Pfaffian to a scalar Toeplitz determinant, 



(34) 



we can 



Pfn) 



pf(T n [a P ]) 

n(ri — 1) 



T„[s. 
T n [s + - 1] 



where I used the fact that TJs. 



1)— P f| 

det(T n [s_]), 

1] = -Tnfa-l^This is the assertion. 



(35) 

□ 



13 0(2n) stands for the orthogonal matrices in E 2 
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3 Subleading order in the NESS EFP asymptotics 

Due to Proposition |8j the study of the large n behavior of the EFP correlation function boils 
down to the analysis of a large truncated Toeplitz operator whose symbol is scalar and has 
the form given in Theorem |6j 

s_(e ifc ) = - (1 -tanh[i(/3-sign(sinA;)5)(A + cos£;)]) , (36) 

see Figure [Q In a true nonequilibrium situation, i.e. for 5 > 0, the r.h.s. of ([36]) is no 
longer continuous. Hence, as described in the introduction, we want to study the asymptotic 
behavior of the EFP NESS with the help of the so-called Fisher-Hartwig theory whose main 



content is summarized in Theorem 26 of Appendix B.3 We first introduce the so-called pure 



jump symbols. For notation and definitions, see Appendix [BTj 

Definition 9 (Pure jump) Let the argument function arg : C \ {0} -> E be defined by z —. 

| z | e iar g ( 2 ) anc / ar g( z ) <= (-• 7r,7r] for all z e C \ {0}. For /3 e C and t e T, the pure jump 
symbol (pp 0jto e PC (T) is defined, for all t e T, by 

m,t (t) := e^ arg (-i) . (37) 



Remark 10 Note that ^ 0)to e PC (T) has at most one jump discontinuity at the point t , 
namely 

m , o (t ±0)=e^. (38) 
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Moreover, the so-called jump phases are defined as follows. 



Definition 11 (Jump phases) Let a e PC (T) with A a = {tj e T\j = l,...,m}, and let 
a(tj ± 0) ^ for j = 1, ...,m. The numbers fa e C for j = l,...,m, called the pure jump 
phases, are defined by 

a{tj - 0) _ ^ 



e 2mp >. (39) 



Next, we define a regularized symbol which will be extracted from s_ below. 

Definition 12 (Regularized symbol) Let tx ■= 1 and t 2 := -1. The regularized symbol 
b P e C(T) is defined by 



Me ») := ( I™Mf J Vii^"'- "°^^> (40) 



rnih) r L (t 2 ) J UsM^), if-ir<k<0, 



r R {ti) 

where, for a = L, R, the function t q : T -> (0, 1) has ffre form 



^(e 1 ") := ^ (1 - tanh [§/3 a (A + cos fc)] ) , (41 ) 



see Figure^ 



Using the pure jump phases and the regularized symbol, we can recast s_ into the fol- 
lowing form. 
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Lemma 13 (Restricted Fisher-Hartwig form) The NESS EFP symbol s_ e L°°(T) has the 

following properties. 

(a) e PC (T) 

(b) A s _={t ll t 2 } 

(c) The jump phases of s at the points t x and t 2 are given by 

A = -^logf^MV (43) 



2vr \T R {t 2 ) j 

(d) The symbol s_ e L°°(T) can be written as 

s_ = bptphfrtpfofr. (44) 

Proof. The assertions (a) and (b) immediately follow from the form of the symbol s e L°°(T) 



given in (36). Moreover, using the choice 



1 _ f s_(tj -0)\ i 



* = *^{r^))-**{mw))' (45) 

where the argument function is given in Definition [9j we get the pure jump phases ((42) and 
fl43~) in assertion (c). As for assertion (d), writing, for k e (-n, n], 

fc+sign( — fc)7r 



^ 1)tl (e i& ) = (^jry ) • (46) 



> = l^J • (47) 

with the sign function defined after ((24), we get equality ((44) in L°°(T) involving the regular- 
ized symbol b P e C(T). □ 



In order to be able to apply the Fisher-Hartwig theory to the symbol s _, we have to make 
sure that the regularized symbol b P e C(T) is indeed sufficiently regular. 

Lemma 14 (Fisher-Hartwig regularity) The regularized symbol b P e C(T) fras f/ie follow- 
ing properties. 

(a) bp(t) ^ fora//t e T 
(<bj ind(6 P ) = 
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(c) bp e Bl(T) 



Proof. It follows from ((40) that b P (t) > for all t e T which implies as serti ons (a) and (b). 
As for assertion (c), we use the fact given in Lemma[27](b) of Appendix p^3] that b P e C l (T). 



This allows us to bound the integrand in d65l) in the Definition 21 of the Besov space Bj(T) 



as 



k 2 H A fcMlLi(T) 



< 



* dk 



d6 / dt \b P (eW +t V) - b P (eW- tk) )\. 

-7T JO 



(48) 



Since, due to Lemma [27] (a)-(c), bp is continuous and differentiable at all but finitely many 
points having a bounded derivative ||&p|| L oo (T) < oo, we have bp e AC(T), and, thus, it follows 
from the fundamental theorem of calculus that b' P e Lip(T) with Lipschitz constant || bp \\ Loo( ^. 
Using the Lipschitz continuity on the r.h.s. of ((45]), we get 



^ \\ A k b p\\ L HT) < 47r2 H & pllL-(T) < °°' 



Hence, we arrive at assertion (c). 



(49) 



□ 



We are now ready to formulate the main result of this note. 

Theorem 15 (NESS EFP asymptotics) For n ->> oo, the NESS EFP behaves asymptoti- 
cally as 



P(n) ~ G(b P ) n n Qp F(s 
where the base of the exponential factor is given by 

'* dk 



Gib, 



(50) 



(51) 



satisfying < G(b P ) < 1 for all inverse temperatures in the range < f3 L < f3 R < oo. 
Furthermore, the exponent of the power law factor has the form 



(52) 



Thus, Q P > if and only if (3 L ^ (3 R . 



Proof. Due to ((44]) of Lemma[T3j the nonequilibrium symbol s_ has the form s = b P (pp utl <pp 2 ,t 2 , 
where t x ^ t 2 and, w.r.t. to the form ([72]) of the restricted Fisher-Hartwig symbol, we have 
«i = «2 = 0, and Re(/3i) = Re(/3 2 ) = from ((42]) and ((43) in Lemma[l3} Hence, assumptions 
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(a) and (b) of Theorem [26] from Appendix [BT3] are satisfied. Moreover, due to Lemma [14} 
assumptions (c)-(d) of Theorem^ are also satisfied by the regularized symbol. Then, ( |50j ) 
follows from ([73} in Theorem [26} and it remains to derive the exponential, the power law, and 
the constant factors in ([50]) with the help of (7ty-(76\. As for G(b P ), we get ([51} from ([74]) 
and ([40]) in Definition [12} Moreover, plugging ([42]) and ((43) into {75]), we find ([52]). Finally, 
using {76]), the last factor on the r.h.s. of ([50]) has the form 



F(s_) = E(b P )2 2 ^ Yl 



bp,+{tj 
bp (tj 



Yl 0(1 + ^)0(1-/3,), 



3=1,2 ^ ' ^ j=l,2 

where E(b P ), b Pt ±(tj), ancl G (! =•= A/) are 9 iv en in (77) , ( [78] ), and ((80]), respectively. 



(53) 

□ 



Remark 16 In Aschbacher (6), I derived a bound on the decay rate of the exponential decay 
for the NESS EFP in the more general anisotropic XY chain. As noted there and discussed in 
the present introduction, Theorem [15] yields that this bound is exact for the special isotropic 
case at handJ3 



A Fermionic quasifree states 

Let 21 be the selfdual CAR algebra from Definition [T} We denote by £(2t) the set of states on 
the C* algebra 210 

Definition 17 (Density) The density of a state u e £(21) is defined to be the operator S e 

£(f)® 2 ) with 0<S* = S<land JSJ = l-S satisfying, for all F, G e f)® 2 , 

u{B*{F)B{G)) = {F,SG). (54) 

A special class of states are the important fermionic quasifree states. 

Definition 18 (Quasifree state) A state u e £(21) is called quasifree if it vanishes on the 
odd polynomials in the generators, and if it is a Pfaffian on the even polynomials in the 
generators, i.e. if, for all Fx, ...,F 2n e f)® 2 and for any n e N, we have 

co(B(F 1 )...B(F 2n ))=pi(n n ), (55) 

where the skew-symmetric matrix fi n e C 2 " x2n is defined, fori,j = 1, 2n, by 

iuj{B(F % )B{F 3 )), ifi<j, 
tt ntij := < 0, if i = j, (56) 

[-u(B(FAB(Ftf, ifi>j. 

14 This can also be seen by directly using Szego's first limit theorem, see for example Bottcher and Silber- 
mann [11 , p.139]. 
15 l.e. the normalized positive linear functionals on 21. 



On the subleading order in the asymptotics of the nonequilibrium EFP 



13 



Here, the Pfaffian pf : C 2 ™ x2n -» C is defined, on all skew-symmetric matrices A e C 2 " x2n := 

{A e C 2nx2n | A t = -A}^\by 

n 

P f ( A ) '■= Si S I1 ( 7r ) II ^T(2j-l),7T(2j), (57) 

7T j=l 

where the sum is running over all pairings of the set {1, 2, 2n}, i.e. over all the (2n)!/(2 n n!) 
permutations n in the permutation group of2n elements which satisfy 7r(2j - 1) < 7r(2j + 1) 
and ir (2 j - 1) < ir(2j), see Figure^ The set of quasifree states is denoted by Q(2l). 

nnn-nmvnrm-rrrifv- 

7T = (123456) 7T = (123546) vr = (123645) vr = (132456) 

Figure 3: Some of the pairings for n = 3. The total number of intersections / relates to the 
signature as sign(7r) = (-l) / . 

The following lemma has been used in Section [2} 
Lemma 19 (Pfaffian) The Pfaffian has the following properties. 

(a) LetX,Y e C 2nx2n withY 1 = -Y. Then, 

pf (XYX*) = det(X) pf (Y). (58) 

(b) LetX g C nxn . Then, 

\ n(n-l) 

J = (-1) ^ det(X). (59) 
Proof. See for example Stembridge p8) . □ 



Pf 








X 




-X 1 






B Toeplitz theory 

The material of this section is taken from Bottcher and Silbermann fTTL fT2j . 
16 A* is the transposition of the matrix A e C nxn . 
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B.1 Function classes 

Let T := {t 6 C | \t\ = 1} stand for the unit circle. We denote by C(T) the continuous 
functions, by C m (T) the m times continuously differentiable functions, and by L P (T) with 
1 < p < oo the usual Lebesgue spaces. Moreover, AC(T) and Lip(T) stand for the abso- 
lutely continuous and the Lipschitz continuous functions on T. Finally, we need the following 
function class. 

Definition 20 (Piecewise continuous) The set of piecewise continuous functions is de- 
fined by 

PC(T) := {a g L°°(T) | The limits lim e ^ + a(e^ k± ^) exist for all k e (-tt, tt]}. (60) 
For a e PC(T) anc/ any t e T of the form t = e ik with k g (— -k, -k], we use the notation 

a(t±0) := lim a(e i(fc±e) ). (61) 

Moreover, the set of jumps ofae PC(T) is defined by 

A a ■= {t g T|a(t-0) ^a(t + 0)}. (62) 
Finally the set of piecewise continuous functions with finitely many jumps is defined by 

PC (T) := {a G PC(T) | card(A a ) < oo}. (63) 



In order to be able to make use of the Fisher-Hartwig theory from Appendix |B.3[ we also 
need to introduce the following function class. 

Definition 21 (Besov space) Let I < p < oo and k e (-7r,7r]. The operator A k : L P (T) ->■ 
L p (T) is defined, on all f e L P (T), and for all 6 e (-tt, tt] , by 

(A fe /)(e i9 ):=/(e^ +fc ))-/(e ie ). (64) 

Moreover, for any n e N, we recursively set A% ■= A fc A^ -1 . For a > and 1 < p < oo, the 
Besov class is defined by 



S«(T) := [/ G L P (T) / dA; |A;|- {1+Qp) ||A™/|| 

J —TT 



V 



Lp(T) 



<oo\, (65) 



where n eN is s.t. n> a 
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Finally, we need the following definition. 

Definition 22 (Index) Let a g C(T) w/Y/7 a(t) ^ for a// 1 g T, and let c : T E w/Y/7 
c g C(T \ {1}) foe s.f. a = |a| e ic . 7/7e index (or winding number) of a is defined by 

. , c(l - 0) - c(l + 0) 

md(a) := ^ '— ^ ^. (66) 

Z7T 



17 Note that the definition does not depend on the choice of such an n. 
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B.2 Toeplitz operators 



For M e N, we denote by £ 2 M (N) the space of all square-summable C M -valued sequences 
Moreover, we set 
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T OO / 



T) := {/ : T — > C MxM | f tj G L°°(T) for all i,j = 1, M}. 



(67) 



f 2 



We then have the following classical result. 

Theorem 23 (Toeplitz) Let {a x } x& c C MxM . The linear operator A : dom (A) c 
(N) defined on all f g dom (A) with maximal domain dom (A) by 

{oo 
i =1 ' i=i 

is a bounded operator on £ 2 M (N) if and only if there exists an a e Lf IxM (T) s.t, forallx g Z, 

/f fro/cfe 

dfc 



— >■ 



(68) 



2tt 



a(e ifc ) e" ifca; . 



Proof. See Bottcher and Silbermann |TTJ p. 186]. 



(69) 

□ 



We then make the following definition. 



Definition 24 (Symbol) Under the assumptions of Theorem\2^ we write the Toeplitz oper- 
ator as T[a]-.= Ae C(£ 2 M (N)). It has the matrix form 



T\a] 



CLq CI— i CL—2 
0,2 CL\ QjQ 



(70) 



The function a e Lf IxM (Y) is called the symbol ofT[o\. If M = 1, the symbol a e L°°(T) = 
L^° xl (T) and the Toeplitz operator T[a] are called scalar, whereas for M > 1 they are called 
block. 

Finally, the Toeplitz operators are truncated as follows. 

Definition 25 (Finite section) Let n e N. The projection P n e C(£ 2 M (N)) is defined, on all 

f := {x u x n , x n+1 , ...} g £ 2 M (N), by 



P n f = {xi, ...,£„, 0,0,...}. 
Moreover, the truncated Toeplitz matrices T n [a] e c MnxMn are defined by 

T n [a] := PnT[a\P n \ TarL (p n ) . 



(71) 



18 W.r.t. the Euclidean norm on 



■>M 
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B.3 Fisher-Hartwig symbols 

The following theorem summarizes the main results on the asymptotic behavior of Toeplitz 
determinants with Fisher-Hartwig symbols. 

Theorem 26 (Fisher-Hartwig) Let a e L°°(T) be a restricted Fisher-Hartwig symbol, i.e., 
for all t g T, the symbol has the form 

m 

a(t) = b(t)H\t-t^ V ^(t), (72) 

and it satisfies the following assumptions: 

(a) ti, ...,t m e T are pairwise distinct points 

(b) aj,(3j e C with |Re(<x,)| < \ and |Re(^-)| < \ for all j = 1, ...,m 
(C) b G L°°(T) n B{(T) 

(d) b(t) ^ for all t g T 

(e) ind(6) = 

Then, forn ->■ oo, the Toeplitz determinant has the asymptotic approximation 

det(T n [a]) ~ G(b) n n Q F(a), (73) 
where the exponential factor and the power law factor are determined by 

G(b) := exp[(log&) ], (74) 

m 

Q ■= E («?-#)■ ( 75 ) 

Here, f x forx e Z denotes the x-th Fourier coefficient of the function f g L^T). Moreover, 
the constant F(a) is given by 

Fin): FdnHbAl,) " 7, ( tj ) - J] , U i 1 ~ t) > < 76 ) 

j=l j=l l<i^j<rn~ i' 
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where we defined 

E(b) := exp[f;/(log6),(log6)_,V (77) 



, i=i 

oo 



b±(tj) := exp X)(log&) ±I <f , (78) 



.i=i 

G °*ft : " 0(1 + 2^) • (79) 

Finally, the function G : C ->• C /s ?r;e enf/re Barnes G -function defined by 

oo 

G(z + 1) := ( 27 r) z / 2 e- 2 ( 2+1 )/ 2 -^ 22 / 2 JJ [(l + -) e"^ 2 /^] , (80) 

n=l 

where 7^ /s Euler's constant. 

Proof. See Bottcher and Silbermann [12, p.582]. □ 



C Regularized symbol 

Since we have to care about the behavior in the neighborhood of the discontinuities, we write 
the derivatives of the regularized symbol explicitly. 

Lemma 27 (Regularity) LetT± := {t e T|Im(t) ^ 0, sign(Im(t)) = ±1}, and set bp J± ■= 



bp\j ± . Then, the regularized symbol b P e C(T) from Definition 12 has the following proper- 
ties. 

(a) 6 PiT± g CJ°°(T ± ) 

(b) b P e C\T) 

(c) The left and right derivatives D±b' p (t) exist for all t e T, but, for j = 1,2, we have 

D -b'p(tj) ^ D + bp(tj). (81) 
Moreover, the second derivative is essentially bounded, 

II & pIIl-(t) < 00 • ( 82 ) 



19 Under the assumptions on b, there exists a logarithm log b e W(T) n b\ /2 {J), where W(T) denotes the 
Wiener algebra, see Bottcher and Silbermann [1 1 , p.123]. Therefore, G(b), E(b), and b±(tj) for j = 1, ...,m are 
well-defined, and the factors on the r.h.s. of are independent of the choice of log b. 
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Proof. Assertion (a) follows from the very form of ((40]). As for assertion (b), we find that the 
one-sided derivatives of b P at the points ti and t 2 coincide, and, for j = 1, 2, are given by the 
expression 



(83) 



Here, for any / : T ->■ C, I used the notation 



D±f{t 1 ) 



1 



lirn l — (f(e^)-f(h)), 

e->0+ ±£ 



D ± f(t 2 ) := lim ±.(f^-))-f(t 2 )). 



(84) 
(85) 



Combining ([83]) with the derivative of 6p,t±, we get 



b' P (e ik ) 



r L (ti)r R (t 2 ) \ 

T R (h)T L (t 2 ) J 



k 

2 IT 



if < k < 



(86) 



where, for a = L, i2, we setf Q (e ifc ) := |(1 + tanh[|/3 a (A + cos &)]). Hence, it follows from ([86]) 
that 6p g C(T). Finally, computing 6 P T± and the left and right derivatives of b' P at the points 
t 1 and t 2 , 



D+b'^h) 
D-b' P (h) 
D + b' P (t 2 ) 
D-b' P (t 2 ) 



\Zr R (t 2 )r L (t 2 ) 

\fT R (t 2 )T L (t 2 ) 



1 . ( r L {t x )r R {t 2 ) \\ 2 , _ A 



2tt H^(*iK(WJ 



/3LT L (t 2 ) 



(87) 
(88) 
(89) 
(90) 



we find, on one hand, that, for j = 1, 2, 



D-^fo) -D+btpfc) = x lT R {tj)r L {tj) [p R f R (tj) - PLfdtj)} 



(91) 



and, on the other hand, we get b P \ f e C(T ± ). Hence, we arrive at assertion (c). 



□ 
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